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Abstract. The Ma-Trudinger-Wang curvature — or cross-curvature 
— is an object arising in the regularity theory of optimal trans- 
portation. If the transportation cost is derived from a Hamiltonian 
action, we show its cross-curvature can be expressed in terms of 
the associated Jacobi fields. Using this expression, we show the 
least action corresponding to a harmonic oscillator has zero cross- 
curvature, and in particular satisfies the necessary and sufficient 
condition (A3w) for the continuity of optimal maps. We go on to 
study gentle perturbations of the free action by a potential, and 
deduce conditions on the potential which guarantee either that the 
corresponding cost satisfies the more restrictive condition (A3s) 
of Ma, Trudinger and Wang, or in some cases has positive cross- 
curvature. In particular, the quartic potential of the anharmonic 
oscillator satisfies (A3s) in the perturbative regime. 



1. Introduction 

Let fi and u be two Borel probability measures on the manifold M 
and let c : M X M — y R be a cost function. The theory of optimal 
transportation begins with the following minimization problem: 

Find a Borel map which minimizes the following total cost among 
all Borel maps ip : M — y M which push fi forward to u: 

(1.1) / c{x, (p{x))d^{x). 

Jm 

Here the push forward ip^fi of a measure yU by a Borel map is the 
measure defined by = for all Borel sets U C M. 

Recently, there have been a series of breakthroughs in understand- 
ing regularity of the solution to this optimal transportation problem. 
After the work of [191 l22l ESI E] , it is clear that regularity of the op- 
timal map is dictated by a quantity introduced by Ma, Trudinger and 
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Wang, called the Ma-Trudinger-Wang (MTW) curvature [12], the cost- 
sectional curvature [TB], or the cross-curvature [H] (see Section [2] for 
the definitions and conventions used in this paper). Efforts have been 
directed at understanding the corresponding condition for smoothness 
and finding examples which satisfy it. Still, little is known about this 
curvature and much of the work has been focused on cost given by the 
square of a Riemannian distance [17] [13] [10] [TT] . 

In this paper, we study the cross-curvature for costs arising from 
natural mechanical systems. More precisely, let (■, ■) be a Riemannian 
metric on the manifold M and let | • | be the corresponding norm. Let 
V : M — y R be a smooth function and consider cost functions c 
defined by minimizing the corresponding mechanical action. Defining 
the Lagrangian v) = — V{x) we let 

(1.2) CTix,y)=M [ Li^{t),j{t))dt, 

Jo 

where the infimum is taken over all smooth curves 7(-) satisfying 7(0) = 
X and 7(1) = y. We recall here that the problem of finding curves 
which achieve the infimum in (II. 2p is called Hamilton's principle of 
least action. We also call the minimizers to the above infimum curves 
of (or paths of) least action. 

In Theorem l3. ll below. we give a characterization of the cross-curvature 
in terms of Jacobi fields. We use it to find a number of examples of costs 
which are non-negatively cross-curved. The first of these is the fam- 
ily of harmonic oscillator actions. In fact, the cross-curvature vanishes 
completely for this new family, which includes the Euclidean distance 
squared cost in case A = 0. 

Theorem 1.1 (Cross-curvature vanishes for the harmonic oscillator). 

Let L : TR" — > R be the Lagrangian defined by 

L{x, v) = 2 1"^!^ ~ 2^ ' 

where A is any symmetric non-positive definite matrix. Then the cross- 
curvature of the corresponding cost defined by /ll.B) is identically zero. 

Non-positive definiteness of the matrix A above is only needed for 
boundedness of the Lagrangian L. One can define a potential V which 
is equal to |x ■ Ax on a large ball B centered at the origin and stay 
bounded outside B. The cross-curvature of the corresponding cost will 
be identically zero in a smaller ball contained in B. 

A posteriori. Theorem 1 1.1 1 can alternately be verified from the explicit 
form of the induced action. For instance when A is negative definite 
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with eigenvalues — A^, — A^, the transportation cost is given by 



Since this action differs from a coordinate reparameterization of the 
standard cost c(x, y) = —x ■ y by null Lagragians — which depend 
separately on x and on y, hence do not affect the minimization (11. ip 
— the vanishing of all cross-curvatures follows immediately from the 
corresponding result for c and the coordinate independence described 
in [16j [25]. However, we originally discovered this vanishing as 
a simple application of Theorem 13. Ij see Section H] below. Absent an 
explicit evaluation of the action integral — which is not possible in 
more complicated settings — it is difficult to guess the cross-curvature 
of a cost function without such tools as our theorem provides. 

The importance of the example from Theorem 1 1.1 1 is the following. In 
the Riemannian case, it is known that non-negative sectional curvature 
is a necessary condition for the cost to satisfy the weak MTW condition 
or to be non- negatively cross curved (see [ISj)- However, the curvature 
term of the Jacobi field equation (13. ip for the above example is given 
by A which is non-positive definite. This shows that the connection 
between the MTW condition and the non-negativity of the sectional 
curvature in the Riemannian case is atypical. 

We also consider the perturbed Lagrangian L^^v) = ||f P — eV{x) 
on the trivial tangent bundle TR". We find computable conditions on 
V for which the corresponding cost defined by (II. 2p satisfies the strong 
MTW condition for all small enough e > 0. More precisely. 

Theorem 1.2 (Perturbed actions which become (A3s)). Assume that 
there exists a constant C > such that 



for all (x, v) in a bounded open subset of the tangent bundle TR" and 
for all unit tangent vectors u, w in the tangent space T^.'RP which are 
orthogonal to each other. Then there is an open set in the product 
R" X R" on which the cost c defined by and satisfies the strong 
MTW condition (A3s) for all small enough £ > 0. 

Theorem 11.21 tells us that the MTW condition is more related to 
change in the curvature along c-segments than to the curvature it- 
self. Note also that the condition in Theorem 11.21 is easily computable 
whereas the cross-curvature is hard to compute in general since — al- 
though it is local in the product manifold M xM [H] — it is nonlocal in 
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M. We further illustrate this computability by considering the radially 
symmetric potentials, i.e. V{x) = f{}^\ In this case, the condition 
in Theorem 11.21 is satisfied if 

(1.3) fi^i+ML) > C> 0, f"{\l±^) > 0, fi^)(\l±^) > 0. 

In particular, the conditions in (11. 3p are satisfied if f{z) = z^, which 
arises from the physical model of the anharmonic oscillator. 

2. Background: Cross-curvature and MTW Conditions 

In this section, we will review some basic facts about the optimal 
transportation problem needed in this paper. The assumptions in the 
theorems stated in this section are simplified to avoid heavy notation. 
The corresponding theorems with relaxed assumptions can be found, 
for instance, in |23j . 

Let M and N be two smooth manifolds (possibly with boundaries) 
and let ^ and v be Borel probability measures on M and A^, respec- 
tively, with compact support. Let c : M x N — y R be a bounded con- 
tinuous cost function, so that the corresponding optimal transportation 
problem is the following: 

Find Borel maps which minimize the following functional among all 
Borel maps (p : M — > N which push /i forward to u (ie. ii{ip~^(U)) = 
u{U) for all Borel sets U C N): 

c(x, ip{x))dn{x). 

The sufficient conditions for the above problem to have a unique 
solution are given by the following theorem (see [23] for the proof). 

Theorem 2.1 (Existence and uniqueness of optimal maps). Suppose 
that the cost c and the measure /i satisfy the following assumptions 

(1) IX is absolutely continuous with respect to the Lebesgue measure, 

(2) the cost c(x, y) is locally lipschitz in x, uniformly in y, 

(3) c is superdifferentiable everywhere, 

(4) the map y i — V dxc{x, y) is injective on its domain of definition. 

Then there is a solution ip ( called the optimal map ) to the above optimal 
transportation problem. Moreover, it is unique fi-almost everywhere. 

In order to discuss regularity of the optimal map, assumptions stronger 
than those stated in Theorem 12.11 are needed. First, we need an open 
set O = AixjV contained in the product MxN for which the following 
assumptions hold: 

(AO): c is C'^ smooth on O, 
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(Al): the map y i — > drcc{x, y) is injective on Af for each x in A4, 
(A2): the map v i — )■ dxdyc{x,y){v) from the tangent space TyN 
to the cotangent space T*M is injective for all (x, y) in O. 
In this paper, we focus on transportation costs given by (11. 2p . The 
above conditions (AO), (Al), (A2) corresponding to these costs are 
discussed in the Appendix. 

Next, we discuss the most important condition in the regularity the- 
ory of optimal maps; called the Ma-Trudinger-Wang (MTW) condition, 
it involves the cross-curvature. To do this, let JCx be the subset of the 
cotangent space T*M defined by 

/C^ = {{x, -dxc{x, y))\y E M} 

and let /C = |J ICx- 

Definition 2.2 (c-exponential) . The c-exponential map exp'^ : /C — > 
M corresponding to the cost function c is defined by 

(2.1) exp'^(x, a) = [dxc{x, ■)]~'^{—a). 

Definition 2.3 (Cross-curvature). The cross-curvature C : TM©/C© 
T*M — > R corresponding to the cost c is defined by 

3 

Cx{u, a, ai) = — 9^9?c(7(t), exp'^(x, a + sai)) 

2 t=s=o 

where 7(-) is any curve which satisfies 7(0) = x, 7(0) = u. Curves of 
the form s 1 — > exp'^(x, a + sai) in the definition of the cross- curvature 
are called c-segments. 



In [Hj, it is shown that the cross- curvature defined above can be 
characterized as the sectional curvature of a certain semi-riemannian 
structure on O in which c-segments form lightlike geodesies — and 
hence the name. The same quantity has also been called the Ma- 
Trudinger-Wang curvature [T2| since it appeared first in Let S be 
the set defined by 

Sx := {(x,M,a,ai) e T^M ® JC^ ® T*M\ai{u) =0}. 

and let 5 = [J S^- In this paper, we use the name Ma-Trudinger- 

Wang curvature to denote the restriction of the cross-curvature C to 
the set S. 

Definition 2.4. The Ma-Trudinger-Wang curvature MTW is defined 
by 

MTW = C 
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Finally, we can state the MTW conditions and the cross-curvature 
conditions. 

(A3w): the cost c satisfies the weak MTW condition on O if 
MTW > on 5, 

(A3s): the cost c satisfies the strong MTW condition on O if it 
satisfies the weak MTW condition on O and MTW^^u, a, ai) = 
only if M = or ai = 0, 

(B3w): the cost c is non-negatively cross-curved on (9 if C > 0, 

(B3s): the cost c is positively cross-curved on O if it is non- 
negative cross-curved on O and Cx{u,a,ai) = only if m = 
or «! = 0. 

The relevance of these conditions to the regularity theory of optimal 
maps can be found in[l9l|16lll71II2|22l|18l[9l[l0l[IIlll2l|8]. Recently, it 
was shown that the conditions (B3w) and (B3s) also have interesting 
applications to microeconomics [7\ and statistics t2]J which have little 
to do with smoothness. 

3. The MTW Curvature and the Jacobi Map 

Let U be an open subset of the cotangent bundle T*M such that all 
elements in U are regular and not on the conjugate locus (the definitions 
of a regular covector and the conjugate locus are analogous to the 
Riemannian case, see appendix for the precise definitions). 

For the rest of the paper, we consider Lagrangians L : TM — t- R 
which are of the form L{x,v) = ^\v\'^ — V{x), where | ■ | denotes the 
norm corresponding to a fixed Riemannian metric. We also identify 
the tangent and the cotangent bundle by this Riemannian metric. 

Let u, V, w be tangent vectors contained in the tangent space T^-M 
at a point x and assume that v is contained in the open set U. Let 
r G R parameterize a family of least action paths t G [0, 1] i — > (p{t, r) 
corresponding to the cost fll.2p with initial velocity v + tw. Let J : = 

dr^p be the corresponding Jacobi field. Let Dt be the covariant 

T = 

derivative along the curve 1 1 — > (p{t, r). Theorem 17. 61 asserts the Jacobi 
field J satisfies the well-known equation 

(3.1) D'^J + R{dt^, J)dt^ + Hess V^{J) = 0. 

All Jacobi fields are solutions to the above second order differential 
equation, the Jacobi equation. We can associate two problems to this 
differential equation: the initial value problem and the boundary value 
problem. We use the term Jacobi map to refer to the map which sends 
the initial value of the boundary value problem to the initial derivative 
(13. 2 p of the corresponding solution whose terminal value is (its first) 
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zero. More precisely, let x and y be two points on the manifold M which 
can be joined by a unique path 7(-) of least action. Let t i — )■ Jit) be the 
Jacobi field along 7(-) such that J(0) = u, J(l) = 0, and J(t) ^ for all 
t in the interval (0, 1). We define the JacoM map J" : MxTM — > TM 
by 



(3.2) 



t=o 



It is not hard to see that the Jacobi map is linear in the variable u. 

The computation of the cross-curvature boils down to the compu- 
tation of the Jacobi map according to the following theorem. In the 
special case of a Riemannian action, Figalli, Rifford and Villani devel- 
oped a related result independently, see Proposition 2.4 of [TT], as we 
learned after the original draft of the present manuscript was complete. 

Theorem 3.1 (Cross-curvature and the Jacobi map). The cross-curvature 
C is given in terms of the Jacobi map J'^ by 



3 

C(m, V, w) = -dl (m, J7''(exp'^(t; + sw),u)) 



s=0 



Here exp'^ denotes the c-exponential (12.11) . 

Proof of Theorem \3.1\ Let t G [0, 1] i — > (f{t,T,s) be a curve of least 
action which starts from the point exp(ru) and ends at the point 
exp^(f + sw). It follows that the cross-curvature (see Definition 12.31) is 
given by 



(3.3) C{u,v,w) 



2 ' ^ 



dt 



=T=0 



Since 1 1 — > (f{t, r, s) is a curve of least action, we have Dtdtf 
as in Theorem 17.61 It follows that 

'1 



dr 



idtifr + viif) 



{Drdtip, dt(p) + dV{dr(p) 

dt{dtip,dr^)+2dV{drip). 



Since the energy i|t>p+\/(a;) is invariant along the curve {x{t),x{t)) = 
{ip, dtf), the left side of the above equation is independent of t. So we 
integrate with respect to t and get 



tdr 



{dr^,dtip) +2 / dV{drV)dt 
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The minimizers t i — > ip{t, r, s) all end at the point exp^(f + sw) 

and the above equation yields 



independent of t. Therefore, drff 
'1 



t=i 



{t - l)dr 



If we set t = 0, then the above equation becomes 



\dM^ + v{^) 



(dr^,dtip) +2 / dV(dr^)dt. 

t=0 



Recall that ip 



exp(ru), so DT-dr^p 



0. Therefore, if we 



i=0 



differentiate the above equation with respect to r, then we have 

9! 



r=0 



{dr^,Drdt^) +2 [ dl{V{^))dt 

T=t=0 Jq 



T=0 



T = t=0 



+ 2 / dfiViip))dt 



r=0 



It follows that the cross-curvature is given by 



Ciu,v,w) = --dX -\dM'-Viip)dt 



2 



S=T=0 

l\dM^ + V{v)-2 [ V{^)dt 
^ Jo 



=T=0 



-dl{u,J\a{s),u)) 



s=0 



□ 



4. A NEW EXAMPLE: THE HARMONIC OSCILLATOR 



In this section, we discuss the example in Theorem 11.11 More pre- 
cisely, we have the following. 

Theorem 4.1 (Cross-curvature of the harmonic oscillator vanishes). 

Let L he the Lagrangian defined by 

Lix, v) = -\v\'^ X ■ Ax 

where A is a symmetric matrix satisfying A < 0. Then the corre- 
sponding cost c defined by ( (i.^) satisfies conditions (AO), (Al), and 
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(A2) on R" X R" (see SectionlEfor the definitions of the conditions). 
The cross-curvature for the cost c is identically zero. In particular, c 
satisfies condition (B3w), a fortiori (A3w), on R" x R". 

Proof. We first show that the corresponding cost given by (11. 2p satisfies 
condition (AO), (Al), and (A2) on R" x R". According to Theorem 
I7.5[ it is enough to check that any two points can be connected by a 
unique minimizers and there are no conjugate points. 

The curves of least action 7(-) corresponding to the Lagrangian 
L{x,v) = ||f P — |x ■ Ax satisfy 

dhit) = -Aj{t) 

by Theorem 17.61 Let ei, ...,e„ be a basis of eigenvectors for the sym- 
metric matrix A. Let x = Yl^=i^^^i — Sr=i'^*^«- 7(0) = ^ 
and 7(0) = f , then 7(-) is given by 



,x ,v )ei 



i=l 



where g^{t, x\v^) = 
Note that the map 



x' cosh(Ait) + sinh(Ait) if ^ 
x' + tv' if Ai = ■ 



v^e'-''{x,v) = J29\^,^' 

i=l 

is a diffeomorphism from the tangent space T^R" to R" for each x. 
Therefore, given any two points, there is a unique path of least action 
joining them. It also follows that there is no conjugate point. 

Next, we show that the cross-curvature is identically zero. The Ja- 
cobi equation for this Lagrangian is given by Theorem 17.61 to be 

dfJ + AJ = 0. 

Note that the matrix A is independent of time t. It follows that the 
Jacobi map J^'^{y, u) depends only on u but not on y. Thus Theorem 
13. II implies the cross-curvature is identically zero. □ 

5. Perturbation by a Gentle Potential 
In this section, we consider the perturbed Lagrangian 

Le{x,v) = ^\v\^ - eV{x) 

defined on the trivial tangent bundle TR", where V : R" — > R is a 
smooth function which is bounded above. We find conditions on V for 
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which the corresponding costs are positively cross-curved or satisfy the 
strong MTW condition for all small enough e > 0. 

The Hamiltonian (after identifying the tangent and cotangent bundle 
using the Euclidean metric) corresponding to the above Lagrangian is 
given by 

He{x,v) = ^\v\^ + eV{x). 
Let e^^^ be the Hamiltonian flow and let $f be the map deflned by 

$f(x,Q;) = TT{e^^'{x,a)). 

Suppose for each e > that Us is an open subset of the tangent 
bundle TR" such that all elements in Ue are regular and outside the 
conjugate locus (see Appendix ^ for the deflnitions). The following 
theorem is a more precise version of Theorem 11.21 

Theorem 5.1 (Gentle potentials yielding positive cross-curvature). 
Let U be a bounded open subset of the tangent bundle TR" such that 



(5.1) / / (M,(l-t)9,'HessV;+i( 

v+sw) 

u) 

Jo Jo 



dtdr > C 

s=0 



for some constant C > 0, for all {x, v) in the set U and for all unit 
tangent vectors u,w in the tangent space T^R". Then, the costs corre- 
sponding to the Lagrangians are positively cross-curved on the set 

K-e = {{x,^l{x,a))\{x,a) E UeDU}. 

for all small enough e > 0. 

If the condition Ii5. 1\) only holds under the assumption that u and w 
are orthogonal, then the costs satisfy the condition (A3s) on the set 
Ke for all small enough e > 0. 

Proof. Let t G [0, 1] i — )■ 'je(t,s) G M be curves of least action corre- 
sponding to the perturbed Lagrangian which satisfy the conditions 
%{0, s) = and dt'je^O, s) = v -\- sw. Let be the Jacobi fleld along 

the minimizer t i — > 7e(t, s) which satisfies Jj 

It follows from Theorem 17.61 that 

dfJs + eRessVyJs = 0. 



u and Jg 



= 0. 

t=i 



■'t "^s 

If we differentiate the above equation with respect to e and let X 



, then we have 

e=0 



deJe , 
e=0 

(5.2) a^X + HessV^gJo = 

with boundary conditions X 



= and X 

t=o 



= 0. 

t=i 
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The family of curves 70 is clearly given by 
7o(t, s) = X + t{v + sw) 
and the field Jq is given by 

Jo(t) = (1 - t)u. 

Therefore, if we integrate the differential equation (15.21) and apply 
the boundary conditions, we get 



(5.3) 



dtX 



t=o 



Jo 



;i -t)HessV; 



x+t{v+sw) 



udrdt. 



Here the integral signs denote componentwise integration. 
Let Ze = Ze{u, V, w) be a smooth function such that 



diJ^^iexp^-^iv + sw),u) 



s=0 



dlJ^°{exY>^°{v + sw), u) + edldtX{u, v, w 



s=0 



+ Zs{u,v,w) 



where lim ^ = uniformly on the bounded set S defined by 

e — p-0 ^ 

S := {(x, u, V, w)\u, v,w E T^M, (x, v) G U, \u\ = \w\ = 1}. 



Since 9f {u, JT'"" (exp^o (f + sw),u)) 



s=0 



dl {u, J'''' {exp"' {v + sw), u)) 



e{u,dldtX{u,v,w)) 



t=s=0 



0, we have 

+ {u,Z,{u,v,w)) 



It follows from ( 15. 3p and the assumptions of the theorem that 



5^ {u, J"' (exp'^-' {v + sw), u)) 



> C" > 



s=0 



for some constant C and for all small enough e > 0. 
Finally, the term 



dl {u, J"' {exp"' (v + sw),u)) 



s=0 



is homogeneous of degree two in both the u and w variables. Therefore, 



dUu,J'^{exp''^{v + sw),u)) > C'\u\^\w\^ 

s=0 

holds on the set 

{(x, u, V, w)\u, v,w E TxM, (x, v) eU}. 

Combining Theorem 13.11 with Definition 12.31 and 12.41 concludes the 
proof. □ 
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6. Radially Symmetric Potentials 

In this section, we simplify the condition in Theorem 15 .ll by assuming 
that the potential V is radially symmetric. More precisely, we have the 
following. 

Theorem 6.1 (Some gentle radial potentials which yield (A3s) costs). 

Setting V = f{\x\'^/2), let U be a bounded open subset of the tangent 
bundle TR" such that 

for all time t in the interval [0, 1], for some constant C > 0, for all 
{x,v) in the setU, and for all unit tangent vectors u,w orthogonal to 
each other in the tangent space T^R"'. Then, the costs corresponding 
to the Lagrangians Li. satisfy the condition (A3s) on the set /C^ for all 
small enough e > 0. 

Proof. Let V{x) = /(^)- A computation shows that the second 
derivatives are given by 

Therefore, the Hessian in this case is given by 

HessV = /'(^)/ + /"(^)x®x, 

where x ^ x denotes the linear transformation defined hj x ® x{y) = 
{x,y) X. 

Assume that \u\ = \w\ = 1 and {u,w) = 0. A computation shows 
that 

v+sw) 

u) 

= + [(a; + tv, twf +{x + tv, tuf]f"'C-^) + 

+ {x + tv, twf {x + tv, uf /(^) (^^) 

>ec. 

Therefore, Theorem 15.11 applies and the result follows. □ 

Remark 6.2. It is clear that the condition /" > C implies that / is not 
bounded above. It follows that there is no function which satisfies the 
conditions in Theorem 16.11 everywhere. 

Example 6.3 (The anharmonic oscillator with a quartic potential well). 
Let / be a function such that /(s) = on [0, K] which stays bounded 
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above on [0,oo[. Then it is clear that / satisfies the conditions in 
Theorem 16.11 on the set 

U = {{x,v) I \x\ + \v\ < K}. 

7. Appendix: Hamiltonian Mechanics on Manifolds 

This appendix is devoted to a discussion of background material 
on Hamiltonian mechanics, and properties of the transportation cost 
functions which arise by minimizing a Lagrangian action, which supply 
some perspective on the results above. Although its contents may be 
familiar to experts, we include it for the readers' convenience, since we 
are not aware of a suitable reference summarizing this material in the 
literature. 

Let L : TM — > R be a smooth function called Lagrangian. We 
define the corresponding cost functions q by 

(7.1) Ci{x,y)=mi L{-f{s),^{s))ds 

Jo 

where the infimum is taken over all smooth curves 7(-) joining x to y 
(i.e. 7(0) = X and 7(t) = y). 

The minimization problem (17. ip which defined the cost functions 
Ct above is called Hamilton's principle of least action. In order to 
characterize its minimizers — called paths of least action — we first 
consider the Legendre transform H : T*M — > R of the Lagrangian L 
defined by 

H{x,a) = sup [a{v) — L{x,v)] 

vGT^M 

where the supremum is taken over all tangent vectors v in the tangent 
space TxM at the point x. The function H is called the Hamiltonian 
corresponding to the Lagrangian L. 

Let g^, be the canonical local coordinates of the cotan- 
gent bundle T*M. The Hamiltonian vector field H of the Hamiltonian 
H is defined via the above local coordinates by 

g_fdH dH _dH _dH\ 
\dpi' ' dpn dq^' ' dq"- J 
The following theorem is classical (see, for instance, |5]). 

Theorem 7.1 (Existence of least action paths). Let (-, ■) be a Rie- 
mannian metric and let \ ■ \ he the corresponding norm. Assume that 
the Lagrangian L satisfies the following conditions: 

(1) the restriction ^ of the Lagrangian L to each tangent space 
TxM has a positive definite Hessian, 
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(2) L satisfies L{x,v) > \v\ + Ci for all tangent vectors v and for 
some constant Ci 

(3) for each compact set K C M , and constant C2 > 0, there is a 
constant C3 such that L{x,v) satisfies L{x,v) > C2\v\ + C3 for 
all tangent vector v in the tangent space T^M and all x in K . 

Then, given any time t > and any pair of points x and y, there exists 
a curve of the form s 1 — > 7(5) := n{e'^^ {x,a)) such that 7(0) = x, 
7(t) = y, and 7(-) achieves the infimum in ( [7. ip . 

Let us consider the optimal transportation problem with cost func- 
tion given by ci defined in (17.11) . Under the assumptions of Theorem 
17. II and that the initial measure fi is absolutely continuous with respect 
to the Lebesgue measure, it is known that there is a unique solution to 
the optimal transportation problem (see [2[[6]). 

Let TT : T*M — y M be the natural projection tt{x, a) = x and let 
$f : T*M — ^ M be the map defined by 

$t(x, a) = Tc{e^"{x, a)). 

Note that the curves of least action take the form t G [0, 1] 1 — > $t(x, a) 
by Theorem 17.11 In the Riemannian case, these extend to geodesies 
which are not necessarily length minimizing and the map $1 is the 
Riemannian exponential map. 

Definition 7.2 (Regular covector). A covector {x,a) in the cotangent 
space T*M is regular if t G [0, 1] 1 — > $t(x, a) is a unique path of least 
action between its endpoints. 

Definition 7.3 (Conjugate locus). A covector (x, a) is in the conjugate 
locus if the map $i(x, ■) does not have full rank at a. 

The conjugate locus can be characterized using Jacobi field as in 
Riemannian geometry. For this, let cr(-) be a curve in the cotangent 

bundle T*M. The vector fields J of the form J{t) = 9,<l>t(a(s)) 

s=0 

defined along the curve 1 1 — > $t(cr(0)) are called Jacobi fields. 

Theorem 7.4 (Characterizing the conjugate locus with Jacobi fields). 

The covector (x, a) is contained in the conjugate locus if and only if 
there is a Jacobi field defined along the curve t \ — > ^t{x,a) which 
vanishes at the endpoints (i.e. J\t=o = = J\t=i). 

Proof. The covector (x, a) is contained in the conjugate locus if and 
only if there is a tangent vector v in the tangent space TaT*M based at 
the point a such that da^i{v) = 0. Let s 1 — > a{s) be a curve in T*M 
such that cr'(O) = v, then $i(a(s)) defines a family of minimizers and 
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J{t) := (9s$t(cr(s)) is a Jacobi field which vanishes at the endpoints. 

Conversely, let J(-) be a Jacobi field which vanishes at the endpoints, 
then there is a curve a{-) in the cotangent bundle such that J{t) = 

ds^ti<^(s)) ■ The vanishing of J(0) implies that cr(-) can be chosen 

s=0 

to be in one cotangent space T*M. The vanishing of J(l) implies that 
the differential of the map $i sends cr'(O) to 0. □ 

T*M 

Let U be an open subset of the cotangent bundle T*M such that 
all elements in U are regular and not on the conjugate locus. Let 
O = Ai X J\f he an open subset of the product M x M contained in 

{(x, <l>i(x, a) e U}. 

Theorem 7.5 (Action minimizing transportation costs). The cost c = 
Ci defined in ( [7. i| ) satisfies the conditions (AO), (Al), and (A2) on 
the set O. 

Proof. We first prove that the function c is smooth on the set O. Let 
(x, a) be a point in the cotangent bundle T*M such that t \ — > $t(x, a) 
is uniquely minimizing and $i(a;,a) = y. Since the differential of the 
map $1 is of full rank on W, there is a neighborhood Ui x U2 in the 
product M xM and a map : U1XU2 ^ T*M such that "^{x, y) = a, 
\E'(a;', y') = x', and $i(a;', "^{x', y')) = y' for all (x', y') in the set Ui x U2. 
By shrinking the set Ui x U2, we can assume that the image "^{Ui x U2) 
is contained in U. It follows that 

c{x',y')= [ Li<l>t{x',^{x',y')),dMx',^{x',y')))dt 
Jo 

for all pairs {x',y') in the set Ui x f/2. It follows immediately that the 
cost c is smooth. 

For the proof of condition (Al), consider the following minimization 
problem: Find a curve which minimizes the following expression among 
all smooth curves 7(-) starting at y 

(7.2) _c(7(l),y)+ / L(7(s),-7(s))rfs. 



Jo 

Clearly, the above functional is non-negative and it is zero if and 
only if s 1 — > 7(1 — s) is curve of least action connecting 7(1) to y for 
(17.11) with Lagrangian L. Let 7(-) be such a minimizer with 7(0) = x. 
Note that the Hamiltonian corresponding to the Lagrangian {x, v) 1 — >■ 
L{x,—v) in (17.21) is given by Hr{x,a) := H{x,—a). Therefore, by [H 
Theorem 2.3] (see also [1]), there exists a covector a in the cotangent 
space T*M such that •jit) = 7i{e^^''{a)) and e^'^''{a) = dxc{'j{l),y). 
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Finally if we let •yit) = 7(1 — t), a(t) = e*^(a;), and ait) = —a{l — t), 
then they satisfy 

7r(e*^(-4c(x,y))) = 7r(e*^(5(0))) = ^{t) 

and 7(-) is a curve minimizing the action (17. ip between its endpoints. 
In particular, we have 

$i(-4c(x,i/))) = 7r(ei-^(-4c(x,y))) = 7(1) 

Therefore, if (x, y) is contained in the set O, then the above equation 
shows that d^c is injective as a map from Af to T*M and (4c) (a) = 

$!(-«). 

The condition (A2) follows from the above characterization of (4c) ~^ 
and inverse function theorem. □ 

In this paper, we focus on Lagrangian arising from natural me- 
chanical systems. More precisely, let (■, ■) be a Riemannian metric 
defined on the manifold M and let | ■ | be the corresponding norm. 
Let V : M — y R be a smooth function on the manifold M called 
a potential. Natural mechanical Lagrangians are Lagrangians of the 
form 

(7.3) L{x,v) = ^\v\^ -V{x). 

In order to apply Theorem 17.1^ we assume that the potential V is 
bounded above. Let 7(-) be a curve of least action corresponding to 
the above Lagrangian (17. 3p and let Df be the covariant derivative along 
the curve 7(-) with respect to the given Riemannian metric (■,■). Let 
R be the Riemannian curvature tensor. Then the minimizer 7(-) and 
the Jacobi fields J(-), in this case, satisfy the following equations. 

Theorem 7.6 (Mechanical geodesies and Jacobi equations). The curves 
of least action 7(-) corresponding to natural mechanical Lagrangian de- 
fined in Iji7.3\ ) satisfy the Newton's second law 

Dtdti = -VV,. 

A Jacobi field J(-) along a least action curve satisfies the following 
equation: 

DlJ + Hess V^{J) + i?(7, J)7 = 0. 
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Proof. Let j^i') be a family of smooth curves such that 7o(-) is a curve 
of least action, %{0) — x, and — y. Then 



• 1 



£=0 



dMt),DA'ye{t) 



£=0 



£=0 



dt 



f 

Jo 



-Dtdao{t)-VV^o{t), deceit) 



6=0 



dt. 



It follows that Dtda = -W^. 

Let s G R parameterize a family of least action curves t 



^{t,s) 



such that ip(t, 0) = j{t). Let J = ds(p . Then, 

= DsDtdtcp - R{dt(p, ds(p)dt(p 
= -DsW^ - R{dt(p, dsip)dtip 
= -Hess V^{ds(p) - R{dt(p, dsip)dt(p. 
After setting s to 0, we get the equation J+Hess Vy{J) + R{'y, J)^ 



as claimed. 



□ 
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